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Abstract 

Convection in a cylindrical container of small aspect ratio is studied. It is known that 
when, in addition to buoyancy forces, thermocapillarity effects are taken into account, 
resonant interactions of two modes may appear. In the case of 1:2 resonance amplitude 
equations are derived, showing the existence of a stable heteroclinic orbit and rotating 
waves, until now not observed experimentally. 
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1 Introduction 



Since the works of Benard |[T|] at the beginning of this century, thermal convection became one of 
the most studied pattern forming systems [0]. In Benard-Marangoni (BM) experiments a fluid 
layer confined in a container open to the atmosphere is heated from below. When the difference 
of temperature between the top and the bottom surface (i.e. AT) is bigger than a critical value 
convection appears. Pearson showed that besides buoyancy forces, thermocapillarity effects 
provide another instability mechanism in this case. 

For shallow layers buoyancy forces dominate, while for thick layers thermocapillarity is the 
main mechanism. Two different regimes appear depending on the liquid depth and the aspect 
ratio a, i.e., the ratio between the horizontal and the vertical dimensions of the vessel. For big 
aspect ratios (a > 20) regular patterns appear at the onset of convection (typically hexagons 
in BM convection). For small a {a < 10) the arising patterns are determined by the geometry 
of the vessel. In a theoretical work Rosenblat et al. studied the linear and weakly nonlinear 
problem with small aspect ratio for cylindrical 0] and square P| containers. In the linear 
case buoyancy and thermocapillary forces were considered. But the extension to the nonlinear 
regime was only made for pure thermocapillary effects. 

There are important experimental contributions in these two regimes. For big aspect ratio 
the most interesting are due to Koschmieder 0, Cerisier and Schatz et al. 0. For small 
aspect ratio Koschmieder and Prahl performed experiments in square and cylindrical con- 
tainers near convective threshold. Ondarcuhu et al. described a series of dynamical patterns 
that look very different in square |T0[ and in cylindrical containers 

We study BM convection in a cylindrical container of small aspect ratio in the usual case 
when both, buoyancy and thermocapillarity forces, act together. In order to study the con- 
vective regime time dependent amplitudes of the unstable modes are considered. Bifurcation 



theory shows |]T2| that in a system with 0(2) symmetry but with the reflection symmetry with 
respect to the middplane broken (like in BM convection or in Rayleigh-Benard convection under 
non-Boussinesq conditions) not only the modulus, but also the relative phases of the amplitudes 
must be taken into account. Codimension two (CT) points, where two modes become simulta- 
neously unstable are specially interesting because they can lead to a resonant interaction with a 
characteristic dynamical behaviour. A similar case has been studied in convection in spherical 



cells (SO (3) symmetry) by Friedrich and Haken |13 
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In Section 2 we present the basic equations and the boundary conditions governing the 
dynamics of the system. In Section 3 we perform a linear analysis and compare our results 
with the experimental ones. Section 4 is devoted to the weakly nonlinear regime. In particular, 
we derive the amplitude equations for a 1:2 resonance, which were also derived for a model 
equation by means of bifurcation theory [l^ and studied in detail in ref. [l^. A physical 



interpretation of these results is made. Finally, we discuss the main conclusions in Section 5. 

2 Evolution equations and boundary conditions 

Consider a fluid layer of depth d confined between a conducting plate heated from below and 
an upper surface open to the atmosphere, under a temperature difference AT. We assume that 
the fluid is Newtonian and that Boussinesq approximation holds. (We will follow the standard 
notation for the parameters of the fluid). 

Then, in nondimensional variables (obtained dividing space, time, velocity and temperature 
by d, d'^/n, K,/d and AT respectively) the equations governing the dynamics of the fluid are: 

V-v = 0. (1) 



dv 

g^+(V.V)V 



Vp + RTe^ + V^v. (2) 



8T 

^+vVT = V2T, (3) 
at 

where the Prandtl number Pr = u/k, gives an estimation of the relative importance of thermal 
and mechanical dissipative processes. In this work we take the limit Pr —>■ oo, equivalent to 
assume that velocity perturbations relax much faster than temperature perturbations. This 
limit allows to neglect both, the advective term and the temporal derivative of the velocity in 
the Navier-Stokes equation, so the velocity field follows adiabatically the variations of tempera- 
ture. The Rayleigh number [R = agd^AT/uK) and the Marangoni number (M = jdAT/ puK,) 
measure the relative importance of buoyancy and thermocapillarity forces. They are not in- 
dependent but related by a constant T = •y/pagd'^ (M = TR). The limit F ^ corresponds 
to Rayleigh-Benard convection while F — oo is reached when gravity forces are absent or the 
thickness of the fluid is very small {R ~ d'^M). In this work we will analyse three situations: 
a) F = 0.05, corresponding to the experiments in refs. |T0|, b) F = 1, for some situations 



analysed in ref. and c) F = 100, when surface tension is much greater that buoyancy effects, 
hke in microgravity. In order to determine completely the problem we must give also some 
boundary conditions (be). The bottom is assumed to be rigid and conducting 

V = 0, T = To at z = 0, (4) 

while for the upper surface, assumed to be nondeformable (which is a good approximation for 
layers which are thick enough) and partially conducting we consider heat transfer and the so 
called Marangoni condition: 

du ^dT dv MdT dT ^ ^ ^ 

— + M— = — + ——=w = — + Bi{T-Tref) = at z = 1, 5 
oz or oz r 0(p oz 

where T^e/ is a reference temperature, (m, f , w) are the components of the velocity in cylindrical 

coordinates (r, 0, z) and Bi stands for the Biot number, a phenomenological parameter that 

accounts for the heat transfer between the fluid and the air. The limit Bi oo corresponds 

to a perfect conducting medium above the fluid while Bi —>■ corresponds to an insulator 

material. We use for the rest of this paper the value Bi = 0.1 typical in experimental situations 
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In a finite container we must consider the horizontal extend by means of another nondi- 
mensional parameter, the aspect ratio a, defined as the ratio between the radius and the height 
of the cylindrical container (a = R/d). To complete the problem we need some bc's on the 
sidewall that depend on the geometry of the vessel. For example, adiabatic walls are described 
by the condition: 

dT 

— — = at r = a. (6) 
or 

In usual experiments [ill plastic sidewalls, matching that be are used. For the velocity, 
instead of the usual non-slip condition we take non-deformable walls on which the tangential 
vorticity vanishes (slippery walls) 

d dw 

u = 7:— (rt>) = — — = at r = a, (7) 
Or Or 

This be is used because the problem become separable, so analytical methods can be used. 
Slippery "walls" are obtained in the separation among convective cells in a big aspect ratio 
system. Although this b.c. does not correspond to that used in experiments, the corresponding 
solutions are expected to reproduce qualitatively the experimental facts. 
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Once the equations governing the dynamics are given we can seek for solutions of these 
equations. In particular, we are interested in those solutions that appear when the trivial one 
(the conductive state) becomes unstable. 



3 Critical modes 



In this section we briefly recall the main results of the linear analysis of the equations and b.c. 
The conductive state has velocity and temperature fields in the form: 



Vcond — 0, Tcond — —Z + Tq 

and conserving only the linear terms we obtain: 

V- v = 0. 

= Pr{-Vn + R9e^ + V^v). 

ot 



de 

di 



w + V^^, 



with be: 



V = = 

du ^Bd dv ^ae de ^ 

— + M— = — + M— = w = — + Bie = 
oz ox oz oy oz 



at 
at 



z = 0. 
z = l. 



(8) 

(9) 
(10) 

(11) 

(12) 
(13) 



de 



d 



dr dr 



dw 
dr 







at r = a. 



(14) 



Solutions of these equations can be written as: v(r,t) = v(r)e'**, e{r,t) = 6'(r)e**, 
7r(r)e**, where s is the growth rate. When marginal stability {s = 0) holds, we have: 

V-v = 0. 

-Vvr + Ree^ + VV = 0. 

\/^e + w = 0. 



7r(r,t) 



(15) 
(16) 
(17) 



Separation of variables allows to consider solutions of the system (|T5|)-(|T^ with be (p^)-(p!^) 
given by: 



Umij[r,(p,z) 



J mi J 



r, 0, z) 



[l/kmi) cos(m0 + (j)m)J'mikmir)DWmij{z) 
-m/kl^-r) sin(m0 + (pm) Jmikmir) DWmij (z) 



[r, 0, z) = cos(m0 + (f)m)Jm{kmir)Wmijiz) 



emijir, 0, Z) = COs(m0 + (l)m)Jm{kmir)Qmijiz) 
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where m = 0,1,2,... is the azimuthal wavenumber, Jm is the Bessel function of order m and 
i is the radial wavenumber, which indexes the values kmi satisfying the condition J'^{ka) = 0. 
Due to symmetry there is no be for the azimuthal equation, so the phase is not fixed and must 
be taken into account. This differs from the analysis of Rosenblat et al [^], in which the phase 
dynamics was not considered. 

The functions Wmij{z), Qmij{z) are solutions of the system 

{D^ - k'^fW - Rk^Q = 0. (19) 
(D^ - k'^)& + W = 0, (20) 

with be 

W = DW = e = at z = 0. (21) 
W = DQ + BiQ = D'^W + Mk^Q = at z = 1, (22) 

where D = -j^. These equations have solutions for several values of R denoted by Rj (M, Bi^ m, i, a) 
so the Rayleigh number can be considered as the eigenvalue of the former system (and j may 
be considered as a vertical wavenumber). The critical Rayleigh number (Re) is the minimum 
eigenvalue for which equations (|T9|)-(^) have a solution. It corresponds to the onset of con- 
vection. Nevertheless we must notice that in a real system we do not fix M but F. Therefore, 
in order to obtain the critical value Rc we eliminate M using the relation M = TR. Then, the 
critical Rayleigh number is: 

Rc = mm Rj{T, Bi, a, m,i) (23) 

mij 

Once Rc determined the rest of the eigenvalues Rj{Mc, Bi, a, m, i) and the eigenfuntions Wmij, 
Qmij are calculated keeping the Marangoni number fixed Mc = TRc. In Fig. 1 we plot Rc 
versus a for the different modes for F = 0.05. For F = 100 we recover the results quoted 
in ref. [Q. Experiments are in qualitative agreement with these results, but a quantitative 
comparison is not possible mainly due to the idealized be considered in our analysis. It was 
noticed that modes with m = 3 appear for a = 5 — 6 in experiments, a factor 2 bigger than our 



calculations ||TT|. Another discrepancy is that the mode m = is observed in convective cells 
with a < 1 a fact that cannot be explained theoretically. However, numerical simulations of 
the Rayleigh-Benard problem with realistic lateral be show the same sequence as in the present 
paper . 
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4 Weakly nonlinear expansion 



We perform a Galerkin-Eckhaus expansion of the fields in terms of eigenf unctions of the hnear 
problem with time-depending amplitudes. After inserting this expansion in the PDE's and 
projecting over the modes of the adjoint problem, an infinite-dimensional system of ODE's for 
those amplitudes is obtained. Then, a center manifold reduction is performed into the few 
modes that remain dynamically active after all transients have relaxed. 

Projecting the nonlinear system on the eigenf unctions of the adjoint system one arrives to 
(see Appendix) 

(l - §) (^^™,^) + (l - ^) «^.^) = ) + (^™.(^ ■ V)^), (24) 

where vj^- • and 0^^, are the adjoint modes and (■ ■ ■) denotes the average over the fluid volume: 



- BdV=— ' ' ' ' ' 
V Jv ira 



{B) = - j B dV = — dz rdr d<p B (25) 



The velocity and temperature fields are expanded in terms of the modes of the linear problem: 

1^1= X^(^mjj (^)e*'"'^ + Amij{t)e ^"^'^)Jm{kmi1^) | q ■'^■^^z) ^ i'^^) 

where Amij = Rmij^-'''^"' is a complex amplitude with an indetermined phase (pm included. We 
restrict our calculations to the neigbourhood of a codimension two (CT) point where two modes 
arise simultaneously. In our case the spatially homogeneous amplitude equations for these two 



modes match the normal form of such a CT point with 0(2) symmetry [12 

Ai=piAi + q[Ar'Al. (27) 

A-m = PmAm + QmAi A^ , (28) 

where Pj,qj are functions of 1^4^^, \Am\'^ and Af^A''^ + A^A''^ satisfying Pj{0) = 0,gj(0) = 0. 
The last terms of these equations are called resonant terms. If m + / — 1 < 3 these resonant 
terms appear in the normal form. 

A strong resonance is obtained between the two modes (11) and (21) for a = acr (see 
Fig. 1), where the value acr slightly depends on F ( acT = 1-15 for F = 0.05, and acT = 1.17 
for F = 1, 100) and the corresponding critical Rayleigh number for this point is denoted as Rct 
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(RcT = 530 for r = 0.05, Rct = 80.5 for T = 1, Rct = 0.894 for T = 100). Under these 
conditions the normal form is: 



Ai = fiiAi + aiAiA2 - aiAil^ip - hiAi\A2\^. (29) 
ia = /U2^2 - Ci2A\ - asAal^aP - 62^2!^!^, (30) 

where we have taken Pi = /ii — ai|v4ip — 6i|v42p + 0{A^), P2 = /^2 — a2|A2p — &2|^iP + 0{A^), 
qi = ai + 0{A^), q2 = -«2 + O(A'). 

The form of Eqs. (p9D-(pO|) does not depend on the details of the system, but arises from the 
resonance and the symmetry properties. The value of the coefficients depends on the specific 
problem and on the truncation order in the Galerkin-Eckhaus expansion. We consider the 
critical modes (11), (21) and a finite set of stable modes with lower growth rate, i.e., ((01), 
(31), (41), (12)), which are then eliminated using the center manifold reduction (see Appendix 
for details). The values of the coefficients under these approximations are given in Table |l|. 

Eqs. (|29|)-(30) have been studied by several authors |14], |15|. Here we summarize the 
main results. Using polar coordinates for the amplitudes Ai = rie**^^ and A2 = r2e*'^^ system 
(|2|)-(|3) is written 

r'l = aiTiTa cos $ + yUiTi — airf — hiVir^. (31) 

r2 = -a2rl cos $ + ^2^2 - a2rl - h2r2rl. (32) 
2 

$ = (^2— - 2air9) sin (33) 
r2 

with $ = 201 — 02- (Rotational invariance allows to remove one of the phase equations). 
There are three different types of steady-state solutions of Eqs. (^)-(|33|): 
1) Pure modes {P±)'- 

n = 0, = — , $+ = 0, = TT. (34) 
02 

These solutions represent convection with only one mode and correspond to a temperature field 
in the form (see Eq. (p6D): 

^(r,0) oc .r?cos(20)J2(A;2ir) (P_) (35) 
V 0.2 

For the sake of experimental comparison we plot the perturbation of the temperature field on 
the upper surface for the pure modes (21) and (11) in Fig. 2, although only the former one 
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is solution of Eq. (PW). The pure mode (21) has two reflection {Z2) symmetries and the (11) 



only the reflection respect to the diagonal that unites the upflow and downflow motions. In 
experiments shadowgraphy and schlieren techniques are used to analyse the patterns |jTO|, [ |TT[| . 
In the corresponding images downflow motions appear as dark zones and upflow motions as 
bright zones. The mode (21) must give images with a cold line (downflow) in a diagonal that 
joins two black zones, while the mode (11) gives up and down motions near the sidewalls on 
the two opposite sites of a diagonal. 

Pure mode solutions exist provided 02/^2 > and lose their stability to mixed modes when 

/ii > ± tti — , (36) 



02 \a2/ 
that we will denote line 11. 
2) Mixed modes (M±): 

= /ii ± air2 - airi^ - 6ir2^, 

= /i2r2 =F asri^ - 62'^iV2 - a2r2^, (37) 

$1 =0, $_ = TT. 



The two mixed modes represent mixed states of modes (11) and (21) (Eq. (p6D ): 

e{r, 0) oc ri cos(0) Ji(A;iir) + r2 cos(20) J2(A;2ir) (M+) (38) 

Two examples are given in Fig. 3. Notice that these two modes mix in such a way that only one 
of the two Z2 symmetries of the (21) mode is broken, i.e, the angle between these two modes 
is (M+) or TT (M_). The mixed mode M_|_ undergoes a Hopf bifurcation on the line given by 
( P7| ) together with 

2 

airi + a2r2 = — , (39) 

2r2 

which gives rise to standing waves (SW) characterized by $ = 0, but time-dependent amplitudes 
ri and r2. We will denote A the line in which this ocurs. 
3) Rotating waves (RW): 

ri = 2— r2 , r2 = , (40) 

«2 A 

_ 1 /i2[2(ai/a2)Qi + ^i] - fii[2{ai/a2)b2 + 02] . 

cos sP — , l^i j 

air2 A 



with 

A = 4— ai + 26i + 2—62 + 02. (42) 
At first order in perturbations they can be represented as: 



^(r,0) (X y^^^^^ \y^2 + + ^o<2(0 + 0it) -<l>)J2(A;2ir)j (/?1^)(43) 

with A, 01 and $ given by Eqs. (0), and (^11) . 

In this state the phases of the two modes 0i and 02 are variable but the relative phase $ is 
constant and different from zero. So, the remaining Z2 symmetry of the (11) mode is broken, 
as it is shown in the two examples given in Fig. 4. These rotating waves exist provided that 
(2/ii + /i2)A > and | cos $1 < 1, i.e., 

2— ai + 61 V2 - ("2— 62 + a2Vil < a?(2/ii + /i2)A, (44) 

OL2 J \ Ot2 /J 

and bifurcate off Mj. at | cos$| = 1 (line fi). 

The pattern of these RW rotates with a phase velocity 



= — 2Q;ir2sin$ = —2a\, 



'\ A V «?A(2/ii+/i2)y'' 

where /3 = 2(a;i/a2)oi + &i and 7 = 02 + 2(Q;i/a2)&2- In Fig- 5 the phase velocity (|45|) is plotted 
as a function of R for several values of the aspect ratio a. This has been calculated from the 
equation 

= (1 - ^^^) , m = 1,2, (46) 

that relates the control parameters /ii, [I2 with the physically relevant nondimensional numbers 
R and a. 

The RW can become unstable to modulated waves (MW) with time dependent amplitudes 
and phase velocities 0i, 02 7^ 0. 

4.1 Heteroclinic orbit 

When condition 

0>^,-':^>-aj!^)^- (47) 



02 \a2 
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is fulfilled, a stable heteroclinic orbit exists. We gathered in Fig. 6a the form of the patterns 
that may appear in four representative points of that orbit. Usually a point of this orbit 
corresponds to a mixture of the modes (21) and (11), but when it reaches the ri-axis the phase 
suddenly jumps in vr. Therefore, the amplitude of the mode (21) reverses sign. When ri ^ 
there is a pure mode (21) that varies its phase in tt/2. Experimentally this heteroclinic orbit 
would lead to patterns evolving from one pure mode to another with a phase changed in 7r/2. 
This orbit is structurally stable to changes in the parameters ai,ai,bi. But it is interesting 
to notice that slight changes in the form of Eqs. (^)-(^) (imperfect symmetry) destroy the 
heteroclinic orbit that degenerates into a limit cycle. In fact numerical noise suffices to convert 
the heteroclinic orbit into a periodic orbit in numerical simulations of Eqs.(^)-(^) shown in 
Fig. 6b. Small thermal noise, cell imperfections, etc, unavoidable in real experiments will 
prevent the formation of the heteroclinic orbit which is replaced by periodic motions between 
modes (11) and (21). 

In Fig. 7 the bifurcation diagram is plotted as function of /ii, fi2 for F = 0.05. This diagram 
can be plotted in the plane {R, a) as shown in Fig. 8. In Fig. 9 the stability diagrams for F = 1 
and F = 100 are represented. The quotations in these figures correspond to the stability regions 
of the different solutions. From Figs. 7 and 8 it can be deduced that for a > acr the dynamics 
begins with a mixed mode (M_) that changes into a rotating wave [RW) when the heating 
is increased. For a < acr convection starts in a pure mode (P) that bifurcates to a mixed 
mode (M_|.). This is replaced by a heteroclinic orbit {H) when the control parameter is rised. 
Notice that {H) is stable in a full region of parameter space. Therefore we hope that a periodic 
dynamics that arises from this heteroclinic connection would be reachable in experiments with 
suitable conditions. 

5 Conclusions 

We performed an analysis of BM convection in a cylindrical container with small aspect ratio. 
For the linear problem we calculated the marginal curve Rc{a) and the unstable modes as a 
function of the aspect ratio a for several values of the ratio between surface tension and buoyancy 
effects (F). These results have been obtained for idealized lateral boundary conditions (slippery 
walls). Therefore, they allow only for qualitative comparison with experiments. The sequence 
of modes and the critical Rayleigh number compares quite well with experimental findings. 



11 



However, the critical aspect ratios calculated are smaller than those observed in experiments. 
This has also been obtained in the case of BM convection in a square container |18|. A possible 
explanation is that the viscous boundary layer in real non-slippery walls rises the aspect ratio 
respect to an idealized slippery one. 

The nonlinear analysis is performed near a codimension two (CT) point where two stationary 
pure modes are simultaneously unstable. We choose a situation where a resonance 2:1 between 
modes (11) and (21) appears. A center manifold reduction is done to determine the coefficients 
of the normal form in the neighbourhood of that point. Depending on the parameter values the 
stationary solutions are: a) a pure mode (P) (with two reflection symmetries), mixed modes 
(M) (with only one Z2 symmetry), rotating waves (RW) (without reflection symmetry) and a 
heteroclinic orbit (H). We calculate the dependence of rotation velocity of RW as a function of 
the order parameter. The heteroclinic orbit is quite robust and appears as a stable solution in 
a big region of the parameter space (see Figs. 8, 9). In an experiment this orbit will degenerate 
into a periodic alternancy between two modes that could be detected in a wide region in the 
parameter space. We hope that new experiments will confirm the presence of a heteroclinic 
conection and of rotating waves in this system. 
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Appendix 



The adjoint linear problem is given by the equations 

VV* - Vvr* + 9*e, = 0. 
V ■ V* = 0. 

V^9* + R*w* = 0, 



(A.l) 
(A.2) 
(A.3) 



with be 



u 
w 

* 

u 



du* 

dz 
d , 



w* = 6 
dv* 



0, ai z = 0, 



(rv 



90* r.-n* n.9w* 

oz oz oz 
dw* 

0, at r = a. 



0, at 



^ ' dr 

This problem is separable and the eigenvalues R^ij 
T5|)-(|r7D. The eigenfunctions are: 



(A.4) 

= 1, (A.5) 

(A.6) 

are the same of those of the system 



<ii = k^iCos{m(t))J'^{kmir)DW:;^i^{z). 
vL,, = -{m/r) sm{'m(t))J^{k^ir)DW:;^i^{z). 
kl 



mi] 



cos{m(j))J^{kmir)W^iAz) 



mi] ^ 



Kij = cos{m(f))Jm{kmir)Q*^ij{z) 
with W^ij and Q^ij solutions of the system 

~|2 K2\a\* I D/^2i 



{D' - k')Q* + Rk'W* 



0. 



{D' - kyw* - Q* 



0, 



(A.7) 



(A.8) 
(A.9) 



with be 



Q* = W* = DW* = 0, at z = 0. 



W* = D^W* = DQ* + BiQ* + k'^M^DW* = 0, at z = 1. 



(A.IO) 
(A.ll) 



These eigenfunctions are introduced to obtain 



de 



1 - ^] «,.^) 



Rr 



'"^^ dt Pr ""'^ dt ' 



4.,(v-v)^) + — (v;;,,.-(vv)v) 



(A. 12) 
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that in the hmit Pr 

R. 



oo leads to 
1 / R 



Rr. 



89. 



V)9). 



(A. 13) 



Introducing expression (|26| ) for the temperature and velocity fields, we obtain an infinite set 
of ODE's for the amplitudes. In order to perform a center manifold reduction we must take into 
account the critical modes and the set of stable modes that appear as cuadratical interaction 
of these ones (so they give thirth order terms in the center manifold equations). This is, the 
set: (Oz), (Ij), (2j), (3z), (4i),z = 1, 2, . . . , j = 2, 3, . . .. Since this is again an infinite set we will 
retain only those modes with the smallest eigenvalues R^ij (see Table 2). Moreover we have 
checked that only modes (01), (31), (41), (12) give an estimable contribution to the coefficients 
of the normal form. With this, the following amplitude equations are obtained: 



Toi^Ol 
TuAu 

T31^31 



QQlllll 1^/^ 



11 + 



0011212. 



12 + 



Q012121 I A 



21 



: eoiAoi + aoioioi(^oi + l^oil ) + 
'- eii^ii + «iiii2i ^11^21^ "111221^12421 + "111101 (^iij4oi + ^11^01) 

+ "111201(^12^01 + ^12^01) + "112131^21^31 + "113141^31^41- 
£^21^21 + "211111 "211212^12 + "211131^11^31 + "211231^12^31 

+ "212101(^21^01 + ^21^01) + "212141^21^41- 

£31^31 + "311121^11^21 + "311221^12^21 + "311141^11^41 + "311241^12^41 
+ "313101(^31^01 + ^31^01)- 

- in Ail + "414101 (A41A01 + A41A01) + "411131^11^31 + "411231^12^31 
+ "412121^21- 

- ei2^12 + "121121^11^21 + "121221^12^21 + "121101 

+ "121201 "122131^21^31 + "123141^31^41- 



with 



> (A.14) 



Ql ■■ 
ijk 



(^mnl^mnl), m = 1 . . . 4, = 1 , 2, 

Rc 
' R 



iWjkl) ( 1 



(6**11^^1)6^, z = 1,2 {t 



R 



■jki 



if ~ (j = 0, 3, 4, A; = 1 or J = 1, A; = 2) 

Oiijklmn = {(^ijli'^kll ' V)6'mnl) 



_ ^ Rc, 

R 

T^Rr 



Rr 



(A.15) 
(A.16) 

(A. 17) 
(A.18) 
(A.19) 



After eliminating the modes (01), (31), (41), (12) with a center manifold reduction, we obtain 
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the equations for the critical modes (11), (21): 



^11^11 — eil^ll + ttiiii2l/i2l/iii — 

I ail2l3lQ;3lil2l _j_ Q111221Q121121 V 4 , , \_/\r,, |2 



121 

^ A £ A _L^, /|2 ( 0212101^012121 _|_ 02121410412121 \ A \ A 

r2l/i2l — 621/121 + 0:211111/1^5^ — J- 1 J- J/l2i|/i2l 



2 f 02121010012121 _|_ 02121410412121 ^ A \ A |2 (A. 20) 



^ 02111310311121 _j_ O2121OIO011111 I O211112O12112I M |2 

fill fn / 21 I 111 • 



^ 631 eol 

In Table the values of these coefficients are given for T = 0.05. Dividing these by Th and r2i 
respectively, we obtain the coefficients given in Table |^. 
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Tables 



r = 0.05 


r = 1 


r = 100 


/ii = .1947ei 


/i2 = .359962 


^1 = 0.06654ei 


^2 = 0.0977162 


fii = 0.0538161 


/i2 = 0.0970862 


Hi = .8376 


2 = .1069 


Hi = 0.1812 


(\2 = 0.01952 


Hi = 0.03681 


02 = 0.002350 


ai = 17.45 


a2 = 6.637 


ai = 2.130 


02 = 1.141 


ai = 0.006209 


a2 = 0.003554 


bi = 10.16 


62 = 28.36 


61 = 1.603 


62 = 4.655 


bi = 0.004855 


62 = 0.01384 



Table 1: Coefficients of the normal form for several values of F. 



Rl^^ = 530.3 


R*^^ = 530.3 


R*ou = 696.4 


i?3]^]^ — 812.2 


Rln = 1326 


R121 — 1332 


i?*2i = 2403 


RI21 = 4039 



Table 2: Eigenvalues of the adjoint problem. 



601 = -3.102 


611 = 1.592611 


621 = 1-473621 


= -.4902 


641 = -1.168 


gi2 = -3.889 


Tn = 8.179 


r2i = 4.093 


"012121 = 7.750 


"011111 = 34.68 


"111101 = -12.77 


"112131 = 9.211 


"112111 = 6.851 


"211111 = -1.666 


"211131 = 4.255 


«2i2ioi = -8.267 


"311121 = -2.726 


"112121 = -1.081 


"212141 = 7.040 


"211112 = -7.594 


f> 111121 = 6.851 


f»iii22i = -11.515 


f» 121121 = 5.685 







Table 3: Coefficients used in center manifold reduction for F = 0.05. 
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Figure Captions. 



Figure 1: Marginal curve representing the critical Rayleigh number Rc as a function of 
the aspect ratio a. The indexes {mi) stand for the azimuthal (m) and radial (i) wavenumbers. 
Only the modes with vertical wavenumber j = 1 are represented because we have assumed the 
ordering Rmn 

Figure 2: Pure modes a) (11) and b) (21). Only the second one appears as a solution of 
Eq. |3^. (Point a in Fig. 8). 

Figure 3: Mixed mode a) M+, b) M_. (Points b and c in Fig. 8). 

Figure 4: Rotating waves for a) $ = 2.0 and b) $ = 2.25. (Points d and e in Fig. 8). 

Figure 5: Phase velocity 02 (in absolute value) as a function of Rayleigh number R when 
a) a = 1.00, b) a = 1.15, c) a = 1.30. The solid line corresponds to the case when the RW 
are stable and the doted line when they are unstable. These RW appear as a supercritical 
bifurcation. From Eqs. (|i5| ) and ( ^6] ) it can be seen that, for e ^ 1, 02 e^- (e = ^^^^^ ) • 

Figure 6: a) Heteroclinic orbit, b) Temporal series of the amplitude ri when the system is 
in the former orbit. After a transient, the heteroclinic conexion merges into a limit cycle due 
to the numerical noise. 

Figure 7: Bifurcation diagram in the plane /ii, fi2- Parentheses indicate unstable solutions. 
The lines Si and S2 represent the limits of existence of the heteroclinic orbit, as given by Eq. 
(|47|). (For further explanations see the text). 

Figure 8: Stability diagrams in the plane {R, a). Only the stable solutions are represented. 
When more than one solution is stable, that with the smallest basin of attraction is represented 
in parenthesis. The points a, b, c, d, e, f correspond to the solutions shown in Figs. 2b, 3a, 3b, 
4a, 4b and 6, respectively. 

Figure 9: Stability diagram in {R, a) space, for a) F = 1 and b) F = 100. The stability 
region for the heteroclinic orbit (H) increses as F increases. 
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Figure 1. 
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Figure 2a. 
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Figure 2b. 
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